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ABSTRACT 

Context. Some T Tauri stars show a pecuhar X-ray spectrum that can be modelled by two components with different 
absorbing column densities. 

Aims. We seek to explain the soft X-ray component in DG Tau, the best studied of these sources, with an outflow 
model, taking observations at other wavelengths into consideration. 

Methods. We constrain the outflow properties through spectral fitting and employ simple semi-analytical formulae to 
describe properties of a shock wave that heats up the X-ray emitting region. 

Results. The X-ray emission is consistent with its arising from the fastest and innermost component of the optically 
detected outflow. Only a small fraction of the total mass loss is required for this X-ray emitting component. Our 
favoured model requires shock velocities between 400 and 500 km s~^. For a density > 10^ cm~'^ all dimensions of the 
shock cooling zone are only a few AU, so even in optical observations this cannot be resolved. 

Conclusions. This X-ray emission mechanism in outflows may also operate in other, less absorbed T Tauri stars, in 
addition to corona and accretion spots. 
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1. Introduction 

Jets and outflows seem to be a natural consequence of disc 
accretion, as they are common phenomena in star forma- 
tion. Jets are an especially active area of study at present. 
Their innermost region is particularly interesting because 
it arises from the deepest in the potential well and thus 
probes the conditions there, which cannot be resolved by 
direct imaging. This innermost region is crucial for our un- 
derstanding of jet launching and coUimation. We need to 
constrain the different theoretical models of stellar winds 



20011) and HH 154 (jFavata et all 120021 : iBallv et all 120031 : 



Favata et al.l |2006| ). X-rays very likely trace the fastest 



Kwan fc Tadcmaru 1988; |Matt & Pudritz 200 ^), X- winds 
Shu et al.lil994) and disc winds (|Blandford fc Pavndll982l: 
[Anderson et al.l f2005) to produce reliable estimates of the 
associated angular momentum loss from the disc and pos- 
sibly the star. This spin-down is needed to explain the slow 
rotation of accreting classical T Tauri stars (CTTS). 

CTTS also have long been k nown as copious emit- 
ters of X-rays (iFei gelson &: Montm erlc 1999), and we now 
recognise that multiple emission mechanisms contribute 
to the observed radiation. Studies of large statistical 
samples clearly show a coronal component, which occa- 



sionally bre aks out in h ard X-ray flares (jPreibisch et al.l 
[2005; Stclz er et al.ll2007[ ). Additionally a hot spot on the 
stellar surface heated by ongoing accretion contributes 
(jGiinther et all I2007D . It shows up in X-ray line ratios 
originating in high-density reg ions, much denser than 



shocks in the jets. The first spatially extended X-ra y emis- 
sion o f a CT TS jet was found using Chandra bv iGiidel et al.l 
()2005l . I2008D in DG Tau. This star shows a spatially re- 
solved, soft emission extending out to 5". It is no t com- 
patible with a point source origin, and lGiidel et al.l (|2008l . 
hereafter GS08) model it by a combination of a linear 
source plus a Gaussian of w 140 A U diameter at the 
maximum distance to the star. Also, iGiidel et al.l (|2007t ) 
find unusual spectral shapes in some CTTS distinguish- 
ing these so-called TAX (two absorber X-ray) sources from 
most other CTTS. Their spectra can be fitted well with a 
highly absorbed hard component and a much less absorbed 
soft component. They suggest that the marginally resolved 
soft excess is emi ssion from the base of the jet (see also 
ISc hneider fc Schm itt 2008). In this paper we develop this 
picture further in order to use the X-rays as a probe of the 
innermost, fastest outflow and to constrain the jet's spatial 
extent and m ass outflow rate. This idea is supported by 
the finding of ISchneider fc SchmittI (|2008f ) that there is a 
small but significant spatial offset between the soft and the 
hard central components in DG Tau, so even the central 
component is marginally resolved. 



SchmittI Hopl 



We focus specifically on DG Tau, since the source is 
studied at many wavelengths, especially the stronger jet 
to the SW, but much less information is available for the 

_ _ _ _ , _ , Argiro ffi et al.l fainter counter-jet to the NE. In the optical the jet has been 

20071) . There is relatively recent evidence of X-rays from imaged with HST/STIS severa l times, revealing not only 
the outflows themselves, starting with HH 2 (jPravdo et al.l jet rotation (jCoffev et al.l 120071 ) . but also the presence of 



any corona (Kastner et al.' '2002; Stclzer 
Schmitt et al. 2005; Giinther et al. 2006; 
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components with differ ent velocities in the a pproaching jet 
with up to 600 km s~^ (jBacciotti et al.|[2000D . after correct- 
ing for projection effects. 

Based on analysis of forbidden lines, it was in- 
ferred that the jet density rapidly falls off with increas- 
ing distance from th e star (jLavallev-Fouquet et al.l [2000t 
iBacciotti et all |2000D The jet extends to 0.5 pc as HH 
702 ( McGroartv et all2007t ). where shock-excited knots are 
clearly seen, but also closer to the star there has to be a 
heating source to sustain the jet emission. Heating by am- 
bient shocks within the material itself is a good candidate. 
iTakami et al.l ()2004l ) found a cool and slow molecular wind 
with an opening angle near 90 deg, wh ich surrounds th e 
fas ter and hott e r atom ic outflow (see also lBeck et al.|[2008l) . 

iRaga et al.l (|200lD presented a 3D hydrodynamical sim- 
ulation of the DG Tau jet. They assume a sinusoidal ve- 
locity variation and precession. The precession broadens 
up the working surfaces in an otherwise well-coUimated 
beam and produces several shocks where the faster ma- 
terial catches up with the previously ejected slower matter. 
In this model the density decreases naturally with the dis- 
tance to the stellar source. 

In this paper we examine the emission measure, col- 
umn density and temperature of the X-ray emitting region 
through spectral fitting (Sect.O and interpret the soft com- 
ponent as a tracer of the fastest component of the outflow 
(Sect. [3]). We then discuss constraints from other wave- 
length bands and their implication for the outflow origin 
in Sect.m We present our conclusions in Sect. [5] and give a 
short summary in Sect. [SI 



2. X-ray observations 



We rely on the X-ray data presented by 'GS08', consisting 
of a 40 ks XMM-Newton observation (ObsID 0203540201) 
and four Chandra exposures (ObsIDs 4487, 6409, 7247 and 
7246), which sum up to 90 ks observation time. We retrieved 
this data from the archives, reduced it using standard SAS 
and CIAO tasks, and analyse it in XSPEC. Because of 
the lo w count rate we merge all available Chandra data. 
iGSOSi show that the error bars on the fitted values over- 
lap between the four exposures for the soft component. 
The luminosity in the hard component varies by a factor 
of three, but the other parameters are compatible. Since 
we are most interested in modelling the soft component, 
variations in the luminosity and temperature of the hard 
component are not important here. We simultaneously fit 
the merged Chandra and the XMM-Newton data with two 
thermal components with individual absorption. We keep 
the parameters for Chandra and XMM-Newton coupled in 
the soft component and leave the temperature and nor- 
malisation between the hard components free, because it is 
observed to be much stronger in the XMM-Newton observa- 
tion due to intrinsic variability. As an illustration we show 
the data recorded by the EPIC/PN camera which covers 
the broadest energy range in Fig. [1] As black and red/grey 
lines we show two fits to the soft components with different 
tempe ratures. The fits shown us e solar abundances accord- 
ing to lGrevesse fc Sauvall (|1998[ ). 

Due to the low count numbers, the statistical uncer- 
tainties are large. In Fig. [1] the black line shows a model 
with kT = 0.3 keV and m oderate absorption, very similar 
to that presented in lGSOSl However, a model with a much 




channel energy (keV) 



Fig. 1. XMM-Newton EPIC/PN data of DG Tau. 
Overplotted are the contributions by the soft component 
for two different scenarios: kT = 0.3 keV ( A^h = 2.3 x 
10^1 cm-2; black hue) and kT = 0.08 keV (A^h = H x 
10^^ cm~^; red/grey line). 



cooler temperature (kT — 0.08 keV), higher emission mea- 
sure and larger absorbing column density reproduces the 
data nearly as well, as illustrated by the red/grey line in 
the figure. Thus, for the soft component, there is an am- 
biguity that makes it difficult to precisely determine the 
temperature, emission measure, and absorption column. 

To quantify this, we have run a grid of models, keeping 
all parameters for the hard component fixed and leaving 
the cool component's volume emission measure (VEM here- 
after), temperature and absorption column density free. 
The resulting ^^-distribution is shown in Fig. [2] At the 
99.95 % confidence level the distribution ranges over ten 
orders of magnitude in emission measure and one order 
of magnitude in temperature, following a narrow valley in 
Temperature- VEM space. The most likely temperature re- 
gion is 0.2-0.4 keV with a second region within the 99% 
confidence contour around 0.08 keV. The secondary peak 
around 0.08 keV is caused mostly by Chandra data; it does 
not appear in a fit to the data from XMM-Newton, which 
is more sensitive to low energy X-rays than Chandra. 

We verified that the position of this valley in x^-space 
does not change if the parameters of the hard component 
are left free. We fit the thermal emission with two APEC 
models. When using MEKAL, the minimum x^-valley is 
in the same place, but the absolute minimum model is 
shifted to lower energies by up to 0.05 keV along the valley. 
We regard this as an estimate for the systematic uncertainty 
in the fit process. 

Due to this and the ambiguity in fitting the temperature 
and VEM of the soft component, rather than taking the 
formal minimum model, we consider a range of param- 
eter space that is statistically allowed by the observations. 
Future observations will be needed to pin down more pre- 
cise values. Thus, we fit a simple power-law to the minimum 
X^-valley of figure [2] This gives 



VEM 



1052c 



0.45 • 



0.33 keV 
kT 



5.5 



(1) 
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Fig. 2. This shows confidence contours (at 68 %, 90 %, 99 % 
and 99.95 % confidence level) fitting the VEM and the tem- 
perature of the soft component. The parameters in the hard 
component are held fixed. Assuming a simple shock model 
(see text for details) the shock velocity is given on the up- 
per horizontal axis and mass loss rates in Mq yr~^ are 
marked in dotted lines in the plot. For some models absorb- 
ing column densities Nn are given in units of 10^^ cm~^. 
The temperatures of the two models shown in Fig. [1] are 
marked with dashed lines. Our fit from Eq. [1] is indicated 
with a red/grey line. 

The emission measure of the minimum x^-valley deviates 
less than a factor of 2 from equation ([T]) in the energy range 
0.1-0.4 keV (Fig. [J). We will use this relationship in the 
following section. 

To get a handle on the physical conditions in the emit- 
ting region further constraints are needed. The fact that the 
central component of the soft e mission does not deviate sig - 
nificantly from a point source (jSchneider fc Schmittll2008l ). 
gives an upper limit on the emitting volume and therefore 
a lower limit on the density. Its size should be less than 1", 
corresponding to 140 AU at the dis t ance of DG Tau. At 
the same time ISchneider fc SchmittI (|2008l ) found the soft 
component to be significantly offset from the hard compo- 
nent at the stellar position by about 30 AU to the SW, the 
direction of the X-ray emitting, approaching jet. 

3. A wind model 

We use a simple, analytical model for the outflow. We dis- 
cuss the observations giving rise to this model in detail in 
Sect. 14.31 Those observations indicate a coUimated inner 
jet, which moves relatively fast. It is surro unded by slower 
and consecutively cooler components (e.g. iBacciotti et al.l 
120001 see Sect. 14.31 for more references). This situation - 
resembling the layers of an onion - is sketched in Fig. [3l 
The outer wind is launched from the disc whereas the ori- 
gin of the inner jets could be either the inner disc region or 
the star itself. Disc winds are launched at a temperatures 
of ~ lO'* K or much less (jArdila et al.l 120021 : iHerczeg et al] 
I2006D . Stellar winds, even if launched at coronal tempera- 
tures, will cool adiabati cally, and th e radiative cooling time 
can also be very short (|Matt fc Pu dritz 2007). Therefore, 
beyond several tens of stellar radii, we expect all compo- 
nents of the outflow to be too cool to emit X-rays. The 
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Fig. 3. Sketch of the geometry. The flow launched from the 
region closest to the star travels with the highest velocity, 
and internal shocks in the flow there can reach the highest 
temperatures. 

outer layer contains molecular hydrogen and radiates pre- 
dominantly in the IR; the more coUimated and faster flows 
car i be observed in th e optical. 

iRaga et aD (|200lf ) modelled successfully the optical ob- 
servations with a jet model using a time-variable outflow 
speed. This gives rise to strong shocks, when faster gas 
catches up with previous ejecta. We assume that in the 
innermost outflow region the velocity of the flow is high 
enough to heat the post-shock material to X-ray emitting 
temperatures. The shock forms sufficiently far from the star 
that the flow is coUimated. 

Using the strong shock conditions we can transform the 
thermal energy kT, where k is Boltzmann's constant and T 
is the temperature, into the pre-shock velocity in the shock 
rest frame Wshock: 



'^shock 



500 km s" 



kT 



0.33 keV 



(2) 



This velocity is given on the upper horizontal axis in Fig.O 
We emphasise that, in the case of a travelling shock wave, 
the pre-shock gas veloc ity is the sti r n of y shnrk and the mo- 
tion of the shock front. iRaga et"aI1 (j2002f ) provided a semi- 
analytical formulation for estimating the cooling length of 
gas heated in fast shocks (Aw ^ 200 km/s). We checked 
this analytic relation with more detailed simulations of 
the s hock cooling zone using our own code (|Giinther et al.l 
I2007D and found that they agree within a factor of a few. 
Thus, for t he semi-ana l ytic t reatment here, we will use the 
formula of iRaga et al.l ( 2002| ). which can be written 



^cool 



20.9AU 



lO^cm 
no 



^shock 



500 km s 



4.5 



(3) 



where no is the pre-shock particle number density, equal to 
a quarter of the post-shock number density n in the strong 
shock approximation. 



4 



H.M. Giinther et al.: Revealing the fastest component of the DG Tau outflow through X-rays 



We assume an idealised cylindrical geometry for the 
post-shock flow, where the length of the cylinder is dcooi- 
The mass flow rate Mghock through the shock front with 
area A is then given by 



M 



shock 



= Av. 



'shockPO) 



(4) 



where po = /iWH^^iono pre-shock mass density of the 

flow, riiono is the initial number density of heavy particles, 
i.e. atoms and ions, /x the mean relative ion particle weight, 
and rriH the mass of hydrogen. Assuming solar abundances, 
we adopt p, = 1.26. It is worth noting that Mghock is less 
than the total mass loss rate in the flow, since tishock is 
measured in the rest frame of the shock, and shocks in 
these systems are observed to travel at speeds of tens to 
hundreds of km s~^. 

The area of the shock can be determined from 



A = ttR^ = 



Volume 



'^cool 

and n\a. 



VEM 

^c^ion'^cool 



(5) 



where Uc and nion are the post-shock electron and ion 
number density. Combing Eqs. [^HSl and assuming rtjon = 



4ni, 



O.SSrie, the mass flow rate can be written 



Mshock « 2.7 • lo- 



ll 



VEM 
1052 cm- 



0.33 keV\ 



1.75 



kr 



(6) 



It is convenient that, by adopting a simple formula for the 
cooling length (Eq. [3]) and a simple shock geometry, the 
determination of the mass loss rate is independent of the 
density of the flow. 

We can now apply this specifically to the soft X-ray com- 
ponent of DG Tau. By combining the observed relationship 
of Eq. [1] with Eq. [SI one obtains 



shock 



1.2 X 10" 



_iiM0 /0.33 keV 



yr 



kT 



7.25 



or, using equation 

^-11 Mq /500 km s 
yr V ^^shock 



Mshock ~ 1.2 X 10" 



-1 



14.5 



(7) 



(8) 



These formulae indicate that the mass loss rate required to 
explain the observed soft X-rays has a very steep depen- 
dence on the observed temperature. It is clear that a better 
determination of the mass loss rate requires a more precise 
temperature measurement. In Fig. [2] lines of constant mass 
loss rate are marked as dotted lines according to Eq. [S] 
Over the full length of the confidence contour the mass loss 
varies by ten orders of magnitude. 

We now look at the physical size of the shock region, R, 
which is the radius of the base of the emitting cylinder of 
shocked gas. From Eqs. ([2]), ([3]) and ([5]), we can solve for 
the radius (in AU), 



R « 0.97 



10^ cm 



/ VEM 



0.33keV\ 



and using the observed relationship for DG Tau (Eq. [T|), 

3.875 

R « 0.65 AU ( ) ( - ) . (10) 



0.33kcVY 



kr 

In Fig. |4] the confidence contours for the lower mass loss 
rates are shown as functions of the radius of the shock at 
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Fig. 4. Confidence contours (at 68%, 90% and 99% con- 
fidence level) are shown for the parameters radius R and 
shock velocity Wshock for the case of a cylindrical cooling 



zone at a fixed density of n = 
the mass loss rates in Mq yr" 



10 cm Dotted lines mark 



the cylinder base and the cooling length assuming a density 
of TT-ion = 10^ cm~^ as motivated from the optical observa- 
tions. For other densites the scalings are R ~ (Eq. ^ 
and dcooi ^ \ (Eq.[3l). 



4. Discussion 

In this section we will use the formulas and condition de- 
rived above to place limits on the physical mechanism pro- 
ducing the observed X-ray radiation. 



4.1. Elemental abundance 

The signal-to-noise ratio is not sufficient to fit elemen- 
tal abundances, so we performed a dditio nal fits with the 
abundances found by IScelsi et"al] (|2007l ) for the Taurus 
molecular cloud. As usual there is a strong degeneracy 
between emission measure and metallicity. On average re- 
duced abundances can be compensated for by a larger emis- 
sion measure. Changes in the composition, that is particu- 
larly the Ne/0 and Ne/Fe ratio, lead to different thermal 
structures. T he confidence cont ours for fits with the abun- 
dances from I Scelsi et al.l ()2007[ ) follow the same valley in 
X^-space, but are displaced about 0.07 keV to higher en- 
ergies. Models below 0.27 keV are excluded on the 99% 
confidence level. 

It is not clear, however, if the coronal abundance is ap- 
plicable to the jet, since the launching region of the X-ray 
emitting, fast jet component is unknown. 



(9) 



4.2. Limits on the size of the shock 

Figure [5] shows two distinct regions of the parameter space 
within the 99% confidence level. That region in the up- 
per left corner corresponds to a scenario of a large emis- 
sion measure, which is heavily absorbed. The tempera- 
ture is only around 0.08 keV and therefore the i^shock ~ 
250 km a~^. As an upper limit we assume a density of 
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10^ cm"^ (see next section). The cooling length accord- 
ing to Eq. [3] should be about 1 AU. In a simple cylindrical 
geometry the volume emission measure, transformed to a 
cylindrical radius with Eq. [9l requires a base shock with a 
radius of 1000 AU, but this would appear as an extended 
source , which is not observed for the central component 
(|GS08l . ISchneider fc SchmittI 12008V Only a density above 
10'' cm~^ could push all shock dimensions below ~ 100 AU, 
so they appear as point source in the X-ray observation. 

We deem this unreasonably high, also the low energy 
region is excluded in fits to the XMM-Newton data alone 
or in models with different - but reasonable - abundances, 
so we concentrate on the less absorbed case with a higher 
Wshock in Fig- m This parameter region can be realised with 
emission from the inner component of the jet only. 

4.3. Limits placed by other observations 

Observations with the HST allow to spatially resolve the 
outflow of DG Tau far below 1" from the central source. The 
forbidden optic al lines observed place density constraints. 
iBacciotti et al.l |2000) conducted seven exposures with a 
long slit parallel to DG Tau's outflow in Ha and forbidden 
optical lines. They analyse the line shift within 0'.'5 from 
the star and find that the faster moving gas is more con- 
fined towards the central axis. The emission of gas faster 
than 200 km s^^ is mostly confined to the inner slit, corre- 
sponding to a cylinder of radius 15 AU or less. Comparison 
of intensity maps in [O i] , [N ii] and [S ii] leads them to 
conclude that the number density in the fast gas is higher 
than 10'' cm~'^ (the line ratios are not sensitive to densi- 
ties above this value) . The lower velocity emission occupies 
a larger area and originates in regions of lower densities 
of the order 10'^ — 10^ cm"-^. So, the inner, dense and fast 
components seem to be surrounded by slower and less dense 
outflows, maybe continuously down to a cool wind observe d 
in H2 with an opening ang le of 90° (|Takami et al.ll2004l ). 
Consistent density estimates of 10^ — 10^ cm"-^ at distances 
less than 0^.^5, deprojected to » 100 AU, were obtained by 
iLavallev-Fouquet et al.l (|2000( ). Strictly speaking, their fine 
ratios give only lower limits on the density for the fastest 
innermost components, too, but in any case the outer layers 
do not contribute significantly to the X-ray radiation. 

From the total line intensity i Lavallev et al.l (jl997[ ) es- 
timate a total mass loss rate of 6.5 • 10~^ Mq yr"-' if the 
gas is heated by shock s. With somewhat similar assumption 
iHartigan et"al] (|1995[) obtain 3 • 10 ''^ yr'^ At a dis- 
tance of 1'.'2 from the central source iLavallev-Fouquet et al.l 
(|2000D estimate a mass loss rate of 1.4- 10"^ Mq yr"^ This 
evidence again favours the lower corner in Fig. [2l shown in 
more detail in Fig. H) 

In optical and IR line emission is detected blue- 

shifted up to deproje cted velo cities of 600 k n i s~^ 

(iLavaUev -Fouauet et al.] I2OOOI : IBacciotti et"al] I2OOOI : 
iPvo et al.. 2003.). An analys is of data covering more than a 
decade bv lPvo et al.l (|2003[ ) indicates that the knots in the 
outfiow move with a common proper motion of 0'.'28 yr~^, 
which translates into a deprojected velocity of 300 km s~^. 

4.4. The shocked gas component 

As shown in Sect. 14.21 the X-ray emission is unlikely to 
originate in the slowly moving components of the out- 



flow. This leaves the innermost component as the probable 
emission origin. Unrealistically high densities are needed to 
produce an emission measure of 10^^ cm~'^ for low shock 
speeds, so we favour a scenario where a strong shock with 
I'shock > 400 km s^^ heats up a relatively small compo- 
nent of the jet. In this case the emission measure is much 
smaller, so the required mass flux to power this emission 
is several orders of magnitude smaller than the total mass 
fiux observed in the optical. The dotted lines in Fig. [4] show 
the mass loss depending on Wshock- Taking into account the 
obs erved proper mo t ion o f the knots, the gas velocities of 
e.g. IBacciotti et al.l (jgOOtf) are too low to explain the X- 
ray spectrum. Because less than 10~^ of the total flow has 
to reach the highest velocity to produce the X-rays, this 
small component was possibly not detected in the optical 
observations. In scenarios with standing (coUimation shock) 
instead of travelling shock waves the innermost component 
is fast enough to explain the observed emission easily. 

After passing through the X-ray shock the material 
quickly radiates away its energy and cools down. Figure ID 
shows all dimensions, the cooling length and the radius of 
the cylinder, to be small, therefore, they are not expected 
to be resolved in the HST observations. The X-ray emitting 
shock is only a small disturbance and within a few AU the 
material has cooled down again to match the temperatures 
of the surrounding gas. The shocks we discuss here cannot 
be resolved in images directly, so only a spectral analysis 
can help to further refine the model. 

4.5. Formation of the shock front 

Although the Chandra observations are separated by nearly 
two years, no e vidence of motion has been detected 
(|Schneider fc Schm itt 2008), but given the size of the er- 
ror circles a velocity Wfront of the shock front of 100 km s~^ 
in the stellar rest frame cannot be excluded. In this case 
the wind mass loss rate M^ind has to be calculated using 
the sum of Wfront and Wshock: 

I'front + ■f^shock .v n 

A/ wind = Aishock (11) 

"^shock 

If these knots represent internal working surfaces the associ- 
ated Wshock is the difference between the gas motion and the 
proper motion about 300 km s^^ for the knots. This by itself 
is not sufficient to power X-ray emission in the high tem- 
perature regime of Fig. [2] The knots seem to be launched 
irregularly, and the later Chandra observations could pos- 
sibly probe a different knot than the earlier observations. A 
relatively slow shock speed of 400 km s~^ corresponding to 
a total mass loss rate of roughly 3 x 10"^" Mq yr~^ (Fig.jl]) 
is then just compatible with the highest optically detected 
outflow velocities. One more scenario can be postulated: 
Instead of a single (stationary or moving) shock front a 
higher number of small internal shocks could produce the 
same X-ray signature. Each of these "shocklets" would not 
show up in the optical observations, because radius and 
cooling lengths were far below the current resolution limit. 
This requires a clumpy outflow, where the jet contains a 
higher number of little knots surrounded by the optically 
detected fast wind. Our calculated shock area and cooling 
volume then represents the sum of all individual shocks. 
This argument is also valid, if the matter passes through 
multiple shock fronts and is reheated several times. The cal- 
culated mass flux rate still represents the summed flow of 
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matter through all shock fronts, but it needs to be divided 
by the number of shocks to obtain the mass loss rate from 
the star. The optical observations discussed in Sect. 14.31 do 
not resolve the postulated X-ray shock front, but the op- 
tically resolved shock fr onts seem to en i erge s eparated by 
150 AU (Al and A2 in iBacciotti et al.l (|2000D ). so it ap- 
pears reasonable that the X-rays may originate in a single 
shock front within 50 AU of the central star. 

From our modelling there is no clear distinction between 
a stationary emission region as provided by a coUimation 
shock or a moving knot. Only further observations with a 
longer time baseline can clarify this point. This is possible 
well within Chandra's lifetime. 

4. 6. The outflow origin 

The X-rays originate in the innermost and fastest com- 
ponent, which may originate fro m the innermost p art o f 
the disc or from the star itself. iMatt fc Pudritd (|2007[ ) 
showed that a T Tauri stellar wind with Mwind ^ 10~^^ 
Mq yr"""^ cannot be cor onally driven, b ut Alfven wave driv- 
ing may be possible (jCranme il I2008D . The spin rate of 
DG Tau makes centrifugal launching of a stellar wind un- 
likely, whereas the outer layers with lower flow velocities 
are mos t likely centrifugally l aunched from the inn er disc 
regions (|Bacciotti et al.ll2002HAnderson et al.ll2003f ). 

4. 7. Relation to extended X-ray jet of DG Tau 

The material heats up while passing through the shock 
front and we could show that the cooling length is very 
small compared to spatial extend of the jet. While cooling, 
the jet looses energy, so the strength of the shock fronts 
decreases, their velocity jump dim inishes with time. In 
their paper Ide Colle &: Ragal (|200()D show this in magneto- 
hydrodynamical models with a time variable ejection veloc- 
ity leading to travelling shock fronts. Additionally, accord- 
ing to the optical observations, the density decreases with 
the distance from the star. Therefore later heating events 
are less luminous, so the resolved source at 5" is weaker. 
Although the resolved shock is a spatially extended source, 
this does not necessarily require multiple heating events 
between 2" and 5", because the coo ling length increases 
dramatically for low densities (jGSOSf ). 

If the cause of the shock is a collision with slower, pre- 
viously ejected matter we can expect this to happen more 
than once. In this picture the knot of X-ray emission at 
5" is caused by matter catching up with an older ejection 
event. [gsosI convincingly show that knots even further out 
are far too faint in X-rays to be detected. 

5. Conclusion 

The soft X-ray emission of DG Tau is consistent with our 
model of a shock in the outflow. The cause of the shock 
is not yet clear, but future observations can reveal if the 
position of the shock front is stationary. In this case it is 
likely caused by the coUimation of the wind into the jet. 
If, on the other hand, the shock front moves, then we can 
take this as indication that the outflow ejection velocity is 
highly variable. In both cases only the fastest components 
of the optically detected outflow provide sufficient energy to 
power the X-ray emission. Possible future grating observa- 



tions could greatly improve our estimate of the temperature 
in the emitting region and narrow down the range of pos- 
sible shock velocities and therefore - via Fig. [4]- the mass 
loss rate. 

DG Tau is an exceptional CTTS in the sense that its 
central component is so much absorbed that the origin of 
the soft emission cannot coincide with the star. This makes 
it an ideal case to search for alternative origins and our 
simple model successfully reproduces the observation. On 
other CTTS, most notably TW Hya, detailed models of 
an accretion shock ex plain the X-rays production very well 
(iGiinther et al.l[200l . It seems possible, that TW Hya also 
has a weak contribution from wind shocks to its emission, 
but its outflows are much weaker than those of DG Tau, 
so the X-rays would be submerged in the stellar emission. 
Next to accretion and coronal activity the wind shocks, as 
discussed in this article, may hold as the third emission 
mechanism of soft X-rays, although it is difflcult to disen- 
tangle the contributions except in special cases. 

6. Summary 

In the X-ray spectrum of DG Tau there is a spatially only 
marginally resolved soft component. Because no grating in- 
formation is available, the temperature is only poorly con- 
strained and plasma models with a wide range of temper- 
atures can reproduce the observations. The absorbing col- 
umn density is negligible in the case of plasma with a ther- 
mal energy of kT = 0.4 keV and rises to N^i > 10^^ cm^^ 
for kT = 0.1 keV. For this second scenario the fitted emis- 
sion measure would be several orders of magnitude higher. 
We reject this solution because it is not stable to small 
changes of the fit parameters and requires physically un- 
reasonable values for the density in the shock region. 

The emission can be explained by a shock and its 
corresponding cooling zone in the innermost and fastest 
component of the optically detected outflow. For densities 
> 10^ cm^'^ this is consistent with all available optical ob- 
servations. Compared to the larger structures seen in the 
forbidden optical lines, this X-ray cooling zone could be a 
relatively thin layer. Only a small fraction of the total mass 
loss is needed in the innermost layer in order to produce the 
observed luminosity. 
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